Resonance Widths in Open Microwave Cavities studied by Harmonic Inversion 
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From the measurement of a reflection spectrum of an open microwave cavity the poles of the 
scattering matrix in the complex plane have been determined. The resonances have been extracted 
by means of the harmonic inversion method. By this it became possible to resolve the resonances in 
a regime where the line widths exceed the mean level spacing up to a factor of 10, a value inaccessible 
in experiments up to now. The obtained experimental distributions of line widths were found to 
be in perfect agreement with predictions from random matrix theory when wall absorption and 
fluctuations caused by couplings to additional channels are considered. 
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Nowadays open wave-chaotic systems are intensely in- 
vestigated experimentally and theoretically (for recent 
reviews see [H, 0, [I])- Especially effects of coupling, ab- 
sorption, and decoherence are of interest. This aspect 
has been studied experimentally in microwave cavities 
0i @i H and acoustics 0], but also in semiconduc- 
tor devices such as quantum dots Till now mainly 
reflection and transmission amplitudes and phases of the 
spectra have been investigated, as it is extremely difficult 
to resolve the resonances in the regime of strong overlap. 
In the case of non-overlapping resonances, i. e. weak cou- 
pling and weak absorption, the width distributions are 
expected to obey a distribution 2] which_has been 
found experimentally in microwave cavities (lO| . There 
is only one experimental work studying resonance widths 
in the overlapping regime by Persson et al. jlli] , where the 
effect of resonance trapping was investigated. In the lat- 
ter work the centered time-delay analysis (CTDA) was 
used to extract the resonance widths. 

Till now an experimental verification of theoretical and 
numerical investigations of the resonance width in de- 
pendence of the number of channels and their coupling 
strengths [13, 13, 14| had been missing. In the regime 



of strong resonance overlap the CTDA method does no 
longer work. In this regime the harmonic inversion (HI) 
method poses an alternative to extract resonances from 
a spectrum. The HI method, introduced by Wall and 



Neuhauser 15( and irnproved by Mandelshtam, Taylor, 
Main and others [l^, 17 1, allows to determine the un- 
known complex frequencies i^n — *7n and amplitudes a„ 
from a time series C(t) — a„ exp[— 27ri(i'„ — ijn)t]- 
The Vn , 7n , and a„ are obtained from a set of nonlinear 
equations which can be solved numerically stable. A de- 
scription of the method for the extraction of resonance 
positions and widths can be found in Ref . [l^. In this Let- 
ter, the method is applied, for the first time, to analyze 
experimental absorption spectra. 

In the experiment a flat microwave cavity of the 
shape of a half Sinai billiard (length a— 43 cm, width 
6=23.7 cm) has been used, with a half disc of diameter 



d = 12 cm attached to the long side. The reflection coef- 
ficient was measured in the range from 1 to 19.4 GHz by 
means of a single wire antenna attached to the cavity. In 
this frequency range the cavity with a height of 7.8 mm 
may be treated as two-dimensional. To improve statis- 
tics the half disc was moved along the wall in steps of 
5 mm to obtain 57 measurements. The complex S matrix 
was measured with a vector network analyzer. Coupling 
strength of the antenna Ta has been determined by the 
mean values of S {Ta = 1 — K-?)]^) and the total absorp- 
tion strength Tt by the decay of the Fourier transformed 
autocorrelation fmiction [19]. For the evaluation of the 
width distribution it showed up to be necessary to split 
Tt into two parts, one constant contribution Tconst and 
a fluctuating one Tfluc, the latter one being mimicked by 
the coupling of Nc fictitious absorbing channels to the bil- 
liard, Tfluc = Nc-Tc 20 . The complete range from weak 
absorption and weak coupling {Tt = 0.75, Ta = 0.116, 4 
to 5 GHz) to strong absorption and nearly perfect cou- 
pling {Tt = 7.0, Ta = 0.989, 14.7 to 15.7 GHz) has been 
covered. The antenna corresponds to a single channel 
since its radius r — 0.2 mm is much smaller than the 
wavelength. The same system had been considered al- 
ready previously in Refs. 0,0], where the reflection co- 
efficient and the Poisson kernel had been investigated. 
The values for Tt in Refs. 0, 0] slightly differ from the 
ones given above, since the fluctuating part of Tt gives 
rise to a deviation from the exponential decay which has 
not been accounted for in the previous works. 

An example of a spectrum in the strongly overlapping 
regime is shown in Fig. [1] as an inset. The spectrum is a 
superposition of Lorentzians, 



^M = i-E 



(1) 



For the application of the HI technique we first have to 
transform the spectrum to the time domain via a Fourier 
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FIG. 1: (Color online) Modulus square \S{t)f of the Fourier 
transform of the part of the spectrum ^(i^) ranging from 14.5 
to 16.5 GHz (see inset). The dashed horizontal line indicates 
the noise level of the Fourier transform. The dotted vertical 
line shows the cut-off time applied for the harmonic inversion. 
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FIG. 2: (Color online) Real and imaginary part of a part of 
the spectrum in the regime of overlapping resonances. Ad- 
ditionally the reconstructed spectrum using the resonances 
identified by the harmonic inversion is shown in dashed line. 
The vertical lines indicate the positions of the resonances and 
the horizontal lines the corresponding widths at half maxi- 
mum. On top of each figure the difference between the mea- 
sured spectrum and the calculated one is shown. The mean 
level spacing A is marked by a horizontal bar. 



transform, 

oo 

S{t) = ^ y e-2--*5(z.)di. = S{tyJ2a^e 



(2) 



for t > 0, and S{t) = for t < 0.^ 

Fig. [D shows an example for obtained from the 

Fourier transform of the part of the spectrum shown in 



the inset. For the HI only the short time domain has 
been used, since the long time tail of the signal does con- 
tain noise only. The cut-off time is marked in Fig. [1] by 
a dotted vertical line. The number N of data points has 
been of the order of a few hundred up to about 1000 
points depending on the frequency range. A matrix of 
rank N/2 is created in the HI procedure 17, Ts*] and its 
eigenvalues and eigenvectors are calculated, yielding N/2 
complex eigenfrequencies Vn — ijn and their residua a„. 
This number usually exceeds by far the real number of 
resonances in this frequency regime, which is approxi- 
mately given by the Weyl formula 



(3) 



where A is the area, S the circumference, ii' is a con- 
stant of order one, and c the speed of light. Weyl's 
formula counts the mean number of resonances below 
the frequency v for closed cavities. In open systems the 
Weyl formula strictly speaking does not apply, but in the 
present situation the error should be small. Wall absorp- 
tion essentially introduces an imaginary part for the res- 
onances, and has only a minor influence on the density of 
states. Only the presence of the antenna may lead to the 
loss of a single resonance due to resonance trapping [ll | . 
To distinguish real from spurious resonances a number 
of criteria have been applied: (i) Resonances close to the 
border frequencies of the Fourier transform window have 
not been taken into account to avoid boundary effects, 
(ii) The HI has been performed twice, first for S{t), and 
additionally for S{t — to) corresponding to a shift of S{t) 
by one data point. Only resonances showing up to be sta- 
ble with respect to this procedure have been taken, (iii) 
Resonances with a width close to the frequency spacing 
dv between neighboring data points have been removed, 
(iv) Resonances with a depth |a„|/7„ exceeding the noise 
level by less than a factor of two have been removed. The 
number of resonances surviving all checks amounted to 
about 80-85% of the value expected from the Weyl for- 
mula. 

In Fig. [2] the real and imaginary part of the evaluated 
spectrum is shown for a small frequency range in the 
regime of overlapping resonances. The solid line corre- 
sponds to the measured spectrum. The extracted res- 
onances are shown as vertical lines and their widths at 
half maximum by horizontal lines. From the extracted 
resonances we reconstructed the spectrum in terms of 
a sum of Lorentzians (see Eq. ([T])). The reconstructed 
spectrum was adjusted by removing an additional linear 
dependence, such that experimental and reconstructed 
spectrum coincide at both boundaries of the window. 
The reconstructed spectrum is plotted with a dashed line 
in Fig. [2] On top the difference between measured and 
reconstructed spectrum is shown, both for the real and 
imaginary part. The main cause for the deviations are 
fluctuations in the experimental baseline. Since the HI is 
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FIG. 3: (Color online) Normalized widths T/A for 14.7 to 
15.7 GHz with Tt = 7.00 and Ta = 0.989. The horizontal 
dashed line corresponds to the width coming from the total 
absorption (Ft = Tt /(47r)), whereas the dotted line (only vis- 
ible in the inset) corresponds to the constant part of the ab- 
sorption (Fconst = T'const/(47r)). The inset is an enlargement 
of the small width region. 



able only to treat superpositions of Lorentzians it cannot 
account for fluctuations of the baseline properly. But the 
local structure is in perfect agreement with the measured 
spectrum. For the frequency range shown in Fig.[2]Weyrs 
formula gives 18.37 resonances. We found 16 resonances 
indicating that we lost two to three resonances, probably 
due to too small amplitudes. Actually the difference be- 
tween original and reconstructed spectrum suggests that 
there are two other resonances at 15.11 and 15.16 GHz. 
They have not been taken into account due to the noise 
cut-off criterion applied. 

To compare our results from an electromagnetic res- 
onator with the quantum-mechanical predictions, we 
mapped the experimental eigenfrequencies Vn and their 
widths 7„ to their quantum-mechanical counterparts by 
En = Vn and F „ — Vnln This procedure is justi- 
fied as long as the linewidths 7 are small compared to 
the frequencies v, which was the case in the experiment 
(7/^ « 10""^). In the effective Hamiltonian approach 
the quantum-mechanical reflection spectrum S{E) can 
be described by 



1 



E-H,, 



(4) 



where W contains the information on the coupling of the 
antenna to the eigenvalues of the closed system ^] . The 
effective Hamiltonian is given by 



(5) 



Width distributions for a random matrix system with 
multiple coupled channels have been calculated by Som- 
mers, Fyodorov and Titov [l4|, reducing for the single 
channel case to 



1 92 /-i 



X^)e-^->^VF{X,y) (6) 
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FIG. 4: (Color online) Normalized width distribution in the 
frequency range 4 to 5 GHz with Ta = 0.116 and Tt = 0.75. 
The black histogram was obtained from the harmonic inver- 
sion, the light one from a Lorentzian fitting procedure of the 
isolated resonances. The dashed curve corresponds to the the- 
oretical prediction from Eq. (O shifted by a constant off-set 
Ft = Tt /(47r) to account for absorption in the walls. The 
solid line in addition takes into account the fluctuations using 
Tfluc = 0.6 and Nc = 10 (see the text for details). The in- 
set shows a comparison between a random matrix simulation 
and the theoretical prediction using the same values for Ta, 
Tt, Tfluc, and Nc as for the experimental data. 
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FIG. 5: Same as Fig. 14] but for the frequency range 14.7 
to 15.7 GHz with Ta = 0.989 and Tt = 7.0. Additionally the 
asymptotics of the theoretical expectation for perfect coupling 
is plotted as a dotted curve (see text). For the simulations 
and the theoretical curve Tauc = 6.4 and A'^c = 20 have been 
used. 



with 



F{X,y) 



ig-X) / dpi - 



dp2 



(A-P2)Vbi-l)(5-P2) 



(7) 



and g = 2/Tq — 1. Here y — — ■£ is negative. In case 
of perfect coupling the tail is algebraically decaying as 
1/(4^). 

Let us now start with the discussion of the experimen- 
tal results on resonance widths. For the part of the spec- 
trum shown in Fig. [2] the resonance width is by far larger 
than the mean level spacing. This can be better seen 
in Fig. [3] where the widths of the resonances are plotted 
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for the frequency range 14.7 to 15.7 GHz in a regime of 
nearly perfect coupling {Ta = 0.989) and large absorp- 
tion {Tf — 7.0) as a function of the mean level number 
calculated from the Weyl formula. All widths have been 
normalized to the mean level spacing A. The horizon- 
tal dashed line corresponds to the contribution to the 
width from the total absorption, Tt/A = Tt/^in) ^ 
(in Ref. 20] there is a factor of 2tt in the denominator, 
resulting from a differing definition of F). A constant 
total absorption just shifts the theoretical curve. This 
has been considered in Figs. [Hand showing the exper- 
imental results in the regimes of weak and strong reso- 
nance overlap, respectively. The dashed lines correspond 
to the theoretical distribution ([7]) shifted by Tt/A. There 
are clear deviations for smaller widths, but the decay of 
the tails are in good agreement with the theoretical de- 
scription, in particular for Fig. [5l The deviations are 
caused by the fact that there is a finite number of ficti- 
tious channels Nc. Instead of using the exact formulas 
of Ref. [14,] for an arbitrary number of channels, which 
lead to serious stability problems in the numerical inte- 
grations, we folded the exact distribution ([7]) with a 
distribution with Nc degrees of freedom and a first mo- 
ment of Ffluc = Jfluc/(47r). The values of Nc and Thuc 
were chosen by minimizing the deviations from the ex- 
perimental distributions. The insets in Figs.|4]and[5]show 
the result of a random matrix numerical simulation, us- 
ing expression ([4]) for the scattering matrix. The same 
number of eigenvalues has been taken as for the experi- 
mental histograms. The solid line corresponds to theory 
using the folding approximation. 

Fig. [4] represents the regime of isolated resonances. 
Here a standard fit procedure in terms of a superposi- 
tion of Lorentzians was still possible. The resulting dis- 
tribution of line widths shown as light line in Fig. [4] is in 
agreement with the result from the HI. A perfect agree- 
ment between experiment and theory is obtained, but it 
is essential to take into account a coupling to additional 
channels. Coupling strength Tq, and wall absorption Tt 
have been obtained independently, as described above. 
Only the splitting of Tt into a constant and a fluctuating 
part has been adjusted, yielding Tfiuc = 0.6, Nc = 10, 

^const = 0.15. 

Fig. [5] shows the result in the regime of nearly perfect 
coupling, where there is a strong resonance overlap, too 
large to allow a direct fit. Here the full potential of the HI 
technique becomes evident. We find a good coincidence 
between experiment and theory in the tail up to a width 
F/A = 12, already without considering the fluctuating 
part of the absorption, but if it is taken into account 
(3fluc — 6.4, Nc = 20, Tconst = 0.6), an agreement over 
three orders of magnitude is found. In the moment the 
physical interpretation of the additional channels is not 
completely clear. If one tries to describe wall absorption 
by coupling a finite number of fictitious channels to top 
and bottom plate, or, alternatively, to the mantle [l9|. 



the resulting channel numbers are by far too large. They 
might be due to capacitances occurring at interfaces be- 
tween the mantle and top and bottom plates. The con- 
stant contribution to the absorption Tconst is of the order 
of magnitude expected from the resistivity of brass. 

In preliminary studies of systems with strong absorp- 
tion we had already been able to analyze resonances up to 
the limit of F/A « 30. This opens a wide variety of pos- 
sibilities to study the properties of the poles of scattering 
matrices in the complex plane 11311. The verification of 
the prediction by Sommers et al. [141] on the distribution 
of resonances is a first illustration in this respect. 
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